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^^ Mode coupling theory (MCT) has had notable successes in addressing the rheology of hard-sphere colloidal glasses, and also 
^N soft colloidal glasses such as star-polymers. Here, we explore the properties of a recently developed MCT-based schematic con- 
^ stitutive equation under idealized experimental protocols involving single and then double step strains. We find strong deviations 
^^ from expectations based on factorable, BKZ-type constitutive models. Specifically, a nonvanishing stress remains long after the 
^H application of two equal and opposite step strains; this residual stress is a signature of plastic deformation. We also discuss the 
QQ distinction between hypothetically instantaneous step strains and fast ramps. These are not generally equivalent in our MCT 
(S| approach, with the latter more likely to capture the physics of experimental 'step' strains. The distinction points to the different 
role played by reversible anelastic, and irreversible plastic rearrangements. 
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1 Introduction 

The nonlinear flow and deformation of a material can reveal a 
wealth of information on its internal relaxation mechanisms. 
The theoretical prediction and understanding of this nonlinear 
rheological response is, however, a formidable task. Ideally, 
one aims to construct a constitutive equation, i.e., a relation 
that predicts the stress tensor at any time ^ as a functional of 
J the deformation historyf^. The task is particularly relevant for 
systems with large internal relaxation times, where the regime 
of nonlinear response is easily reached by experiment. Such 
systems include entangled polymers, in melts and solutions, 
which have long been a prime area of study in nonlinear rheol- 
ogy ^. Recently however, the case of arrested colloidal states, 
comprising either hard or soft spherical particles in suspen- 
sion at high density, has also received considerable exper- 
imental attention^ ^^. These dense colloidal systems have 
relaxation times that are very long (like polymers) or even, 
within the colloidal glass phase, effectively infinite. The non- 
ergodicity of the glass poses special problems to rheological 
theory, and is the cause of important new features such as a 
yield stress. The flow of arrested collo ids ha s recently been ad- 
dressed via phenomenological theoryCSIiil, mesoscopic mod- 
elsl^StZS and first principles methods such as mode-coupling 
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theory (MCT)^^ ^^. An interesting recent development has 
been the formulation of an MCT-inspired schematic consti- 
tutive model capable of addressing arbitrary time dependent 
flows ^^. This is significantly simpler than the full microscopic 
MCT, but appears to capture its main features. In this paper 
we explore some predictions of this schematic MCT approach, 
focusing on single and double step strains. 

Double step strain deformations are a commonly used pro- 
tocol for rheology and materials testing. For simplicity, we 
will only consider shear deformations here, denoting shear 
strain by y(^) and the relevant off-diagonal element of the 
stress tensor by c{t). In such an experiment one first applies a 
(usually large) shear strain yo at, say, time t = 0. Although in 
theory often idealized as a discontinuous step, this will usually 
take some small, finite time 5^ during which a strain ramp is 
applied. After a waiting time At, a second strain yi is applied, 
and the subsequent stress relaxation is measured. Of course, 
the response to a single step already provides information on 
the relaxation processes of the material. However, by tuning 
At to be comparable with the time scale of those processes, the 
double-step strain response can offer a more sensitive probe of 
the dynamics. 

This holds in particular for the case of reversing double step 
strains, where the steps act in opposite directions: yo > but 
yi < 0. The case where the total strain imposed on the sys- 
tem vanishes, yr = yo + yi = 0, allows an especially clear as- 
sessment of the deformation energy and its recovery. Con- 
sider a nonlinear but elastic solid undergoing a single step 
strain: at long times, a residual stress will remain that relates 
to the strain energy stored in the system. Adding a second 
step strain, the residual stresses created by both will combine 
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in a nonlinear fashion. However, for a purely elastic mate- 
rial, or indeed any viscoelastic material with a unique equilib- 
rium configuration (so-called 'anelastic' materials), the resid- 
ual stress must vanish eventually whenever y^ = 0. The resid- 
ual stress will also vanish eventually in any viscoelastic liquid 
(whose longest relaxation time is by definition finite, so all 
stresses relax). 

In view of those comments it is perhaps unsurprising that 
the ultimate vanishing of the stress a{t = ^^^Jt = 0) after 
exactly reversed double-step strain is predicted by most es- 
tablished constitutive models of the nonlinear rheological re- 
sponse in viscoelastic solids and fluids. In particular, the result 
a(oo, 0) = is unavoidable for any constitutive model in which 
the stress may be written as: 



a(0 = ^ y]f{t-t'\yn')ltt'dt' 



(1) 



where y^^/ is the accumulated shear strain between times t 
and t' . This form is (for simple shear flows) equival ent to 
the so-called BKZ (or Kaye-BKZ or K-BKZ) equation'^^^. 
It approximates the response to a general flow by a nonlin- 
ear superposition (called Boltzmann superposition principle^) 
whose kernel is fixed by the stress response \|/(^y) to a sin- 
gle step strain of amplitude y. For double step strain, Eq. ([T]) 
yields 



j(2) 



(0=Cyi(0 
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(t-At) 
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for times after the second step. Here G^^\t) denotes the re- 
sponse after a single step, and G^^\t) the response after two 
steps separated by A^ for ^ > A^ As ^ ^ oo, BKZ-type models 
predict that the response to the double step approaches that of 
a single step of combined magnitude. 

On the other hand, if in a solid material an irreversible plas- 
tic deformations arises between the initially imposed strain 
and its later exact reversal, the residual stress a(oo,0) has no 
reason to vanish. Consider for example a Bingham plastic, a 
hypothetical material of modulus G showing linear elasticity 
up to strain yy and pure plasticity beyond. A step strain of (say) 
3y3;/2 gives a recoverable strain yy; reversing the step takes 
this to -yy/2 so that a(oo) = -Gyy/2. With this in mind, one 
should not be surprised if well founded constitutive equations 
for colloidal glasses, or indeed other materials where plasticity 
plays a major role, end up looking considerably more compli- 
cated than Eq. [T] 

Most constitutive equations to date for viscoelastic liquids 
and solids have been conceived within the framework of con- 
tinuum elasticity theory or fluid mechanics. A notable ex- 
ception that starts from first principles equations of motion is 
the Doi-Edwards theory of polymer melts and solutions, and 
its later refinements^^. Recently, comparable steps towards a 
constitutive model for dense colloidal suspensions have been 



undertaken through an extension of the mode-coupling theory 
of the glass transition base d on an integration-through tran- 
sients formalism (ITT-MCT)'3EIEI1 starting from the Smolu- 
chowski equation including applied flow, and involving a set 
of approximations tailored to capture the collective structural- 
relaxation dynamics of dense liquids, equations are derived 
for the transient nonequilibrium density correlation functions. 
These in turn determine the stress tensor through a suitably 
approximated generalization of a Green-Kubo relation. The 
memory effects incorporated in MCT imply a nonlinear de- 
pendence on the whole strain history. 

ITT-MCT is a microscopic theory. It takes the static struc- 
ture factor of the quiescent system as input, and can there- 
fore predict, in principle, rheological differences between col- 
loidal systems with different interactions (hard and soft col- 
loids; colloids with or without attractive interactions, etc.). 
Unfortunately however, to actually compute its predictions for 
all but the simplest flows remains a technically challenging 
task^^^^^^. To understand better the qualitative and generic 
features of the resulting constitutive equations, a schematic 
model has been proposed that aims to incorporate the essen- 
tial mathematics of the microscopic equations while reducing 
their complexity ^^ . 

In this paper, we analyse this schematic MCT model for 
the case of single and double step strains, focussing in partic- 
ular on the residual stresses following double step strain de- 
formations. We include a careful discussion of how exactly 
step strains should be implemented within the MCT approach, 
bearing in mind that in the laboratory 'step strains' are not 
genuinely discontinuous, but merely rapid on the time scale of 
structural relaxtion. 



2 Schematic MCT 

We now recapitulate the equations defining the schematic 
model for colloidal rheology; for more details and a discus- 
sion of the relations with microscopic ITT-MCT, see Ref.^^. 
For pure shear with a given shear rate y(r), 

a(0= f W)G{t,t'-W)])dt'. (3) 

J — oo 

G{t^t'\ \^{t')]) is a generalized dynamical shear modulus. In 
the nonlinear response regime it depends functionally on the 
strain history y{t' < t) through the accumulated strain eval- 
uated at various times (see below). Outside the steady-state 
regime, it depends (as does any two-time correlation function) 
on times t and t' separately, not just the interval t — t'. In MCT, 
the generalized modulus is expressed through transient density 
correlation functions. For the schematic model, we write 



G(t,t'-[y{t')])=v^^{t,t 



/\2 



(4) 
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Here the parameter Vq sets the scale of the shear modulus; 
for instance the static linear modulus in the quiescent glass is 
Goo = Vq/q where /o is the nonergodicity parameter of the ide- 
alized glass, given by the long-time limit of the density corre- 
lation function (see below). Here we set Va = 100 which gives 
roughly the correct modulus, in units of ksT/d^ with d the 
colloidal diameter, for hard spheres close to the glass transi- 
tionl^^. The density correlator (\>{t^t'), whose dependence on 
strain history is no longer written explicitly but still present, 
obeys a generalization of the Mori-Zwanzig memory equation 



eled schematically sl^^ 



(7) 



where Jff/ = fl, y(x) dx is the accumulated shear strain between 
times t' and t. Here, Jc is a parameter controlling the critical 
strain that is sufficient to break cages. We set y^ = 1/10 to 
match with experimental findings that this deformation is on 
the scale of a typical cage size, given by the Lindemann crite- 



rion' 



,.1Q.27. 



xoM{t,t')^i^{t,t^)^J^ htt'htt"m{t,t")dt"W'.t')dt" = 0. 3 Step Strains within MCT 



(5) 

Here Xq is a microscopic (Brownian) relaxation time; for our 
numerical work we choose time units such that To = 1 . Also, 
m{t^t") is a memory kernel describing the slow dynamics 
close to the glass transition. Following MCT precepts, this 
is approximated schematically as a nonlinear functional of the 
correlator itself. Specifically, following the form adopted in 
the so-called F12 model, we have 



m{t,t') = m[^{t,t')] = vi(^(r,/) + V2(^(^/) 



/\2 



(6) 



Parameters vi and V2 determine the state point of the model, 
controlling in combination the distance from the glass transi- 
tion and its strength (which will depend in principle on the 
nature of the colloidal interactions). In the quiescent sys- 
tem there accordingly appears a line of critical coupling co- 
efficients (v^, V2) that separates liquid-like responses for small 
coupling strengths from solid-like responses at large coupling 
strengths. The solid phase is an idealized glass, in which 
(within MCT) the longest relaxation time is not merely large, 
but infinite. The glass is therefore characterized by a non- 
ergodic contribution to the correlation functions: in a qui- 
escent system possessing an equilibrium Boltzmann distribu- 
tion, ^{t,t') = ^o{t -t') ^ fo>0 for t-t' ^00, The noner- 
godicity parameter /o of the model is the largest non-negative 
solution of 1/(1 — /) = 1 +m[/]. At the glass transition, the 
long-time quiescent limit jumps discontinuously from /o = 
to /o = /^. For our numerical calculations, we follow con- 
ventionE^and set V2 = 2 which implies v\ = 2(>/2— 1) and 
/c = 1 — l/>/2. There remains a single state parameter 8 
measuring how far the sample lies from its glass transition: 
vi = v\ +£/(a/2— 1), V2 = V2. Ergodic liquid states have 
£ < 0, ideal glass states have 8 > 0. 

The "damping functions" h that appear in Eq. ([5]) incor- 
porate the nonlinear reduction of memory effects caused by 
strain. Strain can break the cage around each particle by forc- 
ing rearrangement of its neighbours; this cage is the cause 
of the memory effects which lead to permanent trapping of 
the particle in the quiescent glass state. The damping func- 
tions are decreasing functions of increasing strain, and mod- 



A step strain of amplitude y can, within MCT, be represented 
mathematically in two distinct ways. The obvious method is 
simply to set y(^) = y5(r) in t he eq uations presented above. It 
was the route chosen in Refs.l^^^L, xhis results in an instan- 
taneous, anelastic (or nonlinearly elastic) deformation which 
has no immediate effect on the memory functions: these only 
start to alter once relaxation of the newly applied step defor- 
mation begins. 

The second route is to consider the same total strain to be 
applied at a finite ramp rate. Intriguingly, this is not equiva- 
lent to the delta function even if the ramp rate is later taken to 
be (within the schematic model) infinitely fast. The cause of 
this distinction is subtle, but clear. It stems from the fact that 
in MCT the structural relaxation time in sustained shear is the 
inverse of the shear rate itself. Thus if a large strain y is applied 
at any finite rate, this leads to a plastic deformation whose cu- 
mulative value does not vanish as that rate increases. By the 
end of the ramp, if y is large enough, all memory of the pre- 
vious flow history (including the early part of the ramp itself) 
has been erased, which is not the case in the delta-function 
approach. We will see that consequences of th is for single 
step strains, as considered in previous work^^'^^'are relatively 
minor unless J^Jc- However the additional memory loss 
becomes important in double step strain, where the response 
after a second step depends strongly on how much this erases 
memory of the first. Here we will find qualitative differences 
between the two approaches. 

In deciding which representation is more appropriate, one 
must bear in mind that ITT-MCT, on which the schematic 
model is based, explicitly restricts shear rates to the range of 
small (bare) Peclet number, yxo <C 1. Thus taking the limit 
of a fast finite ramp is, in principle, just as questionable as 
adopting the delta function. However, as will be seen from the 
results that follow, we find no sign of any singular behaviour 
on taking that finite ramp rate to infinity: indeed the results 
are broadly similar to those found for yio — 1 . Thus the fast 
ramp limit may offer a convenient model of ramps that are 
experimentally step-like (fast compared to the observed qui- 
escent structural relaxation times of the material) yet still slow 
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on the microscopic timescale set by Xq. The delta function ap- 
proach, for the reasons given above, does not represent this 
case. It might possibly capture instead ramps where XqJ ^ 1 
(if in consequence plastic rearrangement is not possible); how- 
ever, as stated already, MCT is not to be trusted in that regime. 
Moreover, it is questionable whether conventional rheologi- 
cal experiments on colloidal materials can involve step strains 
that are fast enough to approach this limit of high bare Peclet 
number during the imposition of the step. 

For these various reasons we currently prefer the finite ramp 
approach as a predictive tool for experimental colloid rheol- 
ogy, and develop this next. After that, we summarize for com- 
parison the results of the delta-function approach, before giv- 
ing our conclusions. 



4 ^Step' Strains Comprising Finite Ramps 

Here we consider single and double step- strain experiments in 
which the 'steps' represent strain ramps of finite width 5t. For 
simplicity, we assume a constant shear rate for each of these 
ramps, so that 



y(0 = < 



Yo f or < r < 5^ 

Yi for At^5t <t <At 

else 



-28t 



(8) 



The shear rates are given by Jt = Jt/dt for / = 0, 1. 

In writing the MCT equations for this shear-rate protocol, 
one is led to consider nine different correlators (seven if the 
two shear rates are of equal magnitude) for the calculation of 
G{t). This can be seen from splitting the domain {t' < t) of the 
{t^t') plane into regions bounded by the discontinuities along 
the t- and ^'-axes, and realizing that {t' < 0) is not needed in 
evaluating Eq. ([3]). 

The resulting equations are defered to the Appendix. Let us 
first highlight some important features for the case of a single 
step strain (yi =0). In that case, only three correlators are 
needed. Two of these are well known from previous world^^^ 
the quiescent equilibrium one (|)o(^ — t') (needed for t >t^ > 
80, and the steady-state (|)^^^ {t - 1') (for bt>t> t'). The third 
is a new two-step correlator which we denote ^^^^\t^t')\ this 
is needed for ^ > 5^ > t' . The corresponding G^^^\t^t') enters 
the Green-Kubo integral, 

G^^\t) = f \oG^^^\ty)dt' = yo I G^^^\t,x)dx (9) 
JO Jo 

where we use the overbar on a as a reminder that the ramp 
rate is finite and the superscript (1) to label the single-step 
response. We have also defined G^^^\t^x) = G^^^\t,t') when 



t' =x8t. The equation of motion for the correlator reads 

+/i[yo(l-x)] / h[yo{l-y)]m[^^'''\t,y)]dy(\>,,{{y-x)8t-,%)dy 

+ /z[yo(l-x)] / m\^Q(t-t'')]dtff'^^^^\t\x)dt'' = (10) 
^5/ 



An analysis of the behavior of the solutions to this equation 
at arbitrary 5^ is intricate. Setting aside for a moment the re- 
striction to small Peclet numbers, the limit 5^ ^ is worth 
discussing. Note that the first integral then vanishes as 0{?)t). 
Equation ([T0| thus turns into a one-parameter set of equa- 
tions for X G [0, 1], as stated in the Appendix (Eq. ([25])). As 
explained there, each of these is independently solved by a 



single- step strain correlator (|) 



(1) 

2yo(1-x) 



{t), where ^\q {t) obeys 



xoa4i^(0+4J^(0+^(yo/2)^'m[^o(^-^0]3.'4J^(^0^^'=o. 



k(i)/ 



(11) 



We will meet ^\q {t) again below when we address delta- 
function steps; in fact, it is the correlator found there by setting 
y(0 = To 8(0 1^ ^^^ schematic equations of motion. We thus 
get for ^ > 8^ 



lim a^y}{t) : 
5r^o ^^ 



'I 



yo / dxG, 



2yo(l-x) 



{t-x5t) 



pdx\G^^l{t) (12) 



where Gy^^ {t) is related to the correlator (|)yQ^ (t) via Eq.UI 
The generalization of Eq. ^ to two strain ramps reads 



Jo 



yiG^^^\t-At-ht,t')dt' (13) 



where we have taken the ramps of equal length 8^ and the 
expression holds for ^ > Ar + 28^. Here g\^ ^ is again the shear 
modulus appearing in the single-strain ramp discussed above, 
but evaluated for the strain yi . The equation for G^^^^ is given 
in the Appendix. If one of 8^ or /St remains finite, Eq. ([13]) does 
not reduce to a form that can be expressed through single-step 
response functions. Only if both 8^ ^ and A^ ^ 0, a simple 
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expression can be obtained in analogy to Eq. ([12]), 



bt^o. 






Yi / dxG, 

Jo 

rYo+Yi 



■Yi 



rYO+Yl (]^ n\ (.^ 

= 4^^,^(0-41^(0+41^(^-^0 (14) 

(This is seen by substituting r| = yo ( 1 — x) + yi in the first, and 
r| = yi(l — x) in the second, integral.) In the last equality, we 
recognize the BKZ prediction, Eq. ([2]). It is obtained from the 
ramp-MCT model only in the double limit 5^ ^ and A^ ^ 
(taken in either order), because in this case, cross-correlations 
that appear for finite At become irrelevant. 



4.1 Results for Ramp Approach 

We now discuss numerical results found from Eq. ^ and ([6]) 
using the ramp approach to step strains as detailed above. 
Adapting a standard algorithm for solving the steady- state 
MCT equations, the equations of motion were solved by sepa- 
rately considering the correlation functions defined in the var- 
ious {t^t^) regions between the discontinuities in y (sketched 
in the Appendix). A large time window is covered in the nu- 
merics by repeatedly doubling the initial discretization step 
size. This implicitly assumes that all correlation functions be- 
come slowly varying if the difference of their time arguments 
is large. This is well established for the quiescent and steady- 
state correlators of MCT. Since the equations of motion share 
the same structure, we believe it to hold also for the various 
correlation functions discussed above (and have not yet found 
any exceptions to this numerically). 

We start with the case of a single step strain. In the ideal 
glass phase (8 > 0), the quiescent- state correlation function 
does not decay but exhibits a finite long-time limit /q. As a 
consequence, the memory kernel also does not decay, mo {t -^ 
oo) =mo[fo] = /o/(l — /o) > 0- Moreover a step strain is un- 
able to melt the glass on a permanent basis, so that a finite 
long-time limit exists also for the memory function and result- 
ing correlator (|)(^^) {t^t') needed to compute the stress response 
to such a step, via Eqs. ^ and ([12]). Accordingly, that stress 
response does not decay to zero but to a finite value c(oo). In 
the limit of short ramps, 5^ -^ 0, the result is found as 



lima(^)(oo)=Va / Jodx [ - 

dt^O Jo \ 1 



h[yo{l-x)]mo[fo] 
+ /z[yo(l-x)]mo[/o: 



(15) 




^ approaches an even-function 

- x) as yo ^ oo. The integral is 



For the choice of the damping function in Eq. (|7]), this inte- 
gral can be evaluated; for yo ^ oo it gives lim^f^Qa^^\oo^ -^ 



Fig. 1 a(oo) after step strain with amplitude yo, for the F12 model 
at the glass transition, 8 = 0, following a constant- shear-rate ramp of 
width 5^ (dash-dotted: analytical result for 5^ -^ 0; circles: for 3t as 
labeled; squares: for fixed yo = yo/S? as labeled). The result a(oo) 
predicted using a delta-function step strain is also shown (dashed 
curve). 



Va(7i/4)yc/o > 0, i.e., the large-yo behavior saturates to a con- 
stant. The dash-dotted curve in Fig. [T] shows the analytical 
result of Eq. ( p3] ) together with numerical results for finite 5^ 
(solid lines with circle symbols). In each case, a linear regime 
extends up to yo ~ yc, indicating the response of an elastic 
solid, a = Gooyo, where Goo is the glass plateau modulus. Af- 
ter that, a sublinear regime indicates plastic deformation. 

The large-yo asymptote can be understood by recognizing 
that in Eq. ([H]), yo/i[yo(l -^ 
representation of (7r/2)yc5(l 
therefore dominated by x = 1. In other words, as 5^ -^ 0, the 
shear rate % ^ ^, and the shear-induced loss of memory is 
perfect for all t' <bt. Hence the the t > dt response is un- 
affected by the imposed shear except for a narrow window at 
the time t^ = 5^ where the shear rate is just being switched off. 
This is in contrast to the 8-step model, explored (for a single 
step) in previous full^and schematic ^^ ^^ MCT calculations, 
whose result for the residual stress a(oo) is also shown in FigIT] 
(dashed line). Here, the limit yo ^ oo yields zero, approached 
as oc yo/z(yo/2)^ ^ I/Tq 1^ ^^^ ^^^e. 

As previously discussed, the finite-ramp approach becomes 
questionable when yo -^ oo, since the microscopic derivation 
of ITT-MCT assumed small Peclet numbers, yoXo <^ 1- How- 
ever, applying the same strain in ramps with lower shear rates 
shows that the fast ramp limit is nonsingular; this is shown 
by the green squares in Fig. [T] For these slower ramps, after 
an elastic regime for yo ^ yc, a constant asymptote is again 
quickly approached for large yo. This is because the dynamics 
in the interval [0, dt] is governed by the transient steady-state 
correlator, which decays on a time scale Xy^ ~ Jc/jo and be- 
comes a function only of yo = yo^ only at long times. For 
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Fig. 2 Infinite-time stress plateau a(oo) after reversing double step 
strains in the glass (F12 model, 8 = 10""*), as a function of the first 
and second step sizes, yo positive and yi negative. Top: Time 
between the strains A^ = 10, with strain ramps of width 5^ = 1. 
Bottom: At = 6t = 0A. 



Yo ^ Yc, the ramp time is always long compared to this relax- 
ation, 5^ :^ XjQ . The integral in Eq. ^ then extends over the 
full steady-state correlator and no longer scales with yo, lead- 
ing to a Yo -independent constant. In consequence, for Yo ^ Yc. 
the residual stress a(^)(oo) is the same as would arise after 
cessation of steady shear. (Note that, within our schematic 
MCT, the latter quantity attains a finite limit as the shear rate 
diverges.) 

We now turn to the case of a double step strain, where the 
residual stress is due to a nonlinear combination of the two 
single steps. Fig. [2] shows results for the resulting residual 
stress a^^) (oo) after a reversing step strains, where Yo > but 
Yi < 0, as a function of the two strain amplitudes. 

Figure [2] shows a remarkable difference of the MCT model 
to the expectation from BKZ-type constitutive equations when 
At is large enough. In the latter, the nonlinear superposition 
principle, Eq. ([T]), implies that the response to an arbitrary 
strain history can be constructed from those to single step- 
strain experiments. In particular, Eq. ([2]) implies a^^\^) = 
GyQ:^y^(oo), so that the plot in Fig. 2 should be inflection- 
symmetric around the line yo = — Yi . As a consequence, ap- 
plying two large strains in opposite directions but such that the 




Fig. 3 Stress recovery after exact step strain reversal: a(oo) as a 
function of and normalized by the elastic response Gooyo for waiting 
times A^ = 0.01, 0.1, 1, 100, and 1000 (top to bottom). F12 model, 
8 = 10~^. Solid lines: ramp model with bt = I. Dashed lines: a(oo) 
for 5- step model. 



total strain y > 0, results in a positive residual stress according 
to BKZ-type models, no matter how large the waiting time be- 
tween the strains. In the ramp-step MCT model, a physically 
more plausible prediction emerges. Whereas for A^ ^ and 
5^ -^ 0, the nonlinear superposition principle is formally re- 
covered, cf. Eq. ( p4| ), for A^ ::^ To. the residual stress a(oo) is 
not a function of yo + Yi alone. Instead it is seen from Fig. [2] 
that the result becomes independent of yo if |Yi| ^ Yc- The 
physics was already described above for a single step: when 
the strain amplitude of the second step is large, \yi\dt ^ Yc 
holds. The system then approaches a steady- state flow dur- 
ing the second strain ramp, which erases all memory of the 
past deformation history, including everything related to the 
first ramp (no matter how large was yo). Thus, we obtain for 
Yi > jc the result a}^]y^{oo) -^ ay|^(oo). Only if the second 
step is small can there be any positive residual stress from the 
first step. The BKZ limit of Eq. ([14]) is only reached if both 
A^, 5^ <C To- even for arbirarily high shear rate, Eq. ^ implies 
that To governs the fastest relaxation possible for 4>(^^0- ^^ 
both the time between the ramps and the duration of the sec- 
ond ramp are short compared to this relaxation, some memory 
can be kept of the first step. 

We next specialize to the case of exactly canceling double 
step strains, yi = — yo. In Fig. [3] we show the residual stress 
normalized by the linear elastic stress, a(oo)/(GooYo) as a func- 
tion of Yo, using a ramp time 5^ = 1 chosen such that the bare 
Peclet number is modest (< 1) within the parameter range ad- 
dressed. The results contrast strongly the class of constitutive 
equations based on the BKZ form Eq. (|2]), for which all the 
curves in Fig. |3] would be identically zero. In the linear re- 
sponse regime, yo -^ 0, the ratio a(oo)/yo does not vanish but 
attains a limiting value that depends on A^ Around yo ~ Yc. 
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Fig. 4 Time-dependent stress relaxation a{t) after a double strain 
ramp, with bt = 0.01, F12 model for 8 = 10~^. The first strain has a 
magnitude yo = 0.1, and the second yi = —Jo/2 (label B), 
respectively yi = — yo (label C). Solid lines: MCT predictions. 
Dashed lines: BKZ prediction, using MCT's single- strain ramp 
result as an input to Eq. ([2]). 



a minimum occurs that becomes more pronounced if At is in- 
creased, until it saturates for time delays beyond the relax- 
ation time of the single step- strain response. At large strains 
Yo > Yc, a(oo) approaches the same saturating value as follow- 
ing a single step, which shows up here as a slowly decaying 
1/Yo asymptote at large strains that is independent of Al (The 
delta-step results, shown in Fig.|3]as dashed lines, will be dis- 
cussed later.) 

Figure |4] shows the time-dependent stress relaxation after 
reversed step strains, for the two cases Yi = — Yq /2 (s ome- 
times classified as type B in the rheology liter atureEED), and 
Yi = — Yo (type C). A strain ramp with 5^ = 0.01 was cho- 
sen, and an initial strain of Yo = 0.1. We note that G{t) can 
be a non-monotonous function, seen here for the type B curve 
which exhbits a maximum at ^ ^ 1 . Such nonmonotonic vari- 
ation arises when the relaxations following the two steps have 
different time scales. According to the discussion of Fig. [2j 
such an interference of the two strain ramps is only possible if 
At and 5^ are comparable to Xq- 

Dashed lines in Fig. [4] show the BKZ prediction accord- 
ing to Eq. (|2]). The BKZ factorization overestimates a{t) 
in type B and similar partially reversed double- step strains. 
This is a known effect in polymer melt rheology ^^^^, where 
qualitatively similar results have been found and the devia- 
tion has been attributed to yisco-anelastic effects ^^ or irre- 
versible breaking of networks (i.e., a persistent loss of mem- 
ory)^. In our model, the difference arises because the system 
keeps memory, but an analysis is not so simple due to the com- 
plicated structure of our constitutive equation. The deviations 
from BKZ are not particularly marked for small and interme- 



diate times. At large times, the BKZ approximation for type C 
differs qualitatively from our model, as explained in connec- 
tion with Fig. [3] Note that in a liquid, this difference can only 
be seen if the quiescent relaxation time is large enough. 



5 Step Strains as Discontinuities 

We now discuss how the ramp-step model discussed so far 
differs from a treatment of instantaneous step strains as delta 
functions within the schematic MCT approach. 

For a double step strain, one can formally set y(0 = yo8(0 + 
yi8{t-At). The total strain Sitt > At is hence y = Yo + Yi • In 
formulating the equations for this idealization, there appear 
damping functions of the form hto = h[yto] with 



Jto = Yo 



Jo 



x)(ix + Yi fovt>At 



(16) 



The integral in this equation can take on any value between 
zero and unity, depending on the chosen sequence of functions 
used to construct the distribution 5(x) in the limit 5^ ^ 0. The 
fact that Eq. ([16]) is formally undefined is of course itself a 
warning that step strains in MCT are not what they seem. 

One obvious choice is to include the full strain in the inte- 
gral, hence y^o = Yo + Ti • Another common interpretation is to 
set the integral to 1/2 (interpreting it as a unit step function 
whose value at the jump discontinuity is set to the midpoint in 
reference to Dirichlet's theorem of Fourier representations). 
We will make the latter choice in what follows (but mention 
the differences where these are marked). The midpoint inter- 
pretation was already used in previous MCT calculations for 
single step strains ^^ ^^, although there the factor 1/2 can any- 
way be absorbed by a redefinition of y^. 

In Eq. ([3]) our choice leads to 






;yi(o 



■-yoG^^Ut)- 



-yi4!^(^- 



-At) 



(17) 



valid for t > At. The second term is the response to a single 
delta-step of magnitude yi , as defined previously via Eq. ^ 
in terms of the correlator 4>yo (0 th^t obeys Eq. (U) . The first 
term G^^\t) cannot be expressed in terms of single-step re- 
sponse functions alone and therefore, as with the ramp-step 
approach, we are dealing with a constitutive model that is 
not of BKZ factorizable form. The corresponding correlator 
^YoJi (0 obeys Eq. ([26]) evaluated at x = 1/2 (see Appendix). 
Before describing the numerical results based on these 
equations we consider first the limit At ^ 0. Then, the first 
integral in Eq. ([26]) can be neglected, and that equation attains 
the form of Eq. ( [TT] ). We get 



lim a^ 

M^O 



,(2) 
Yo,Yi 



(0 



^yoC2Yi(0- 



-yi4!^(^- 



-At) 



(18) 
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Note that this is not equivalent to a single step of the com- 
bined magnitude. In contrast to the BKZ form, Eq. ([2]), the first 
term is not obviously split into two contributions of the form 
G^^\t). Note also how the combination yo/2 + yi enters in the 
first term. This is a direct result of our midpoint treatment of 
the integral in Eq. [16] With this choice only, it is easily es- 
tablished that for an exactly reversed step strain (yi = — yo), as 
At -^ 0, the stress response vanishes identically for all strain 
amplitudes. (Interpreting the integral of Eq. ([16]) instead to 
yield yo + yi , this cancellation would only occur inside the 
linear-response regime.) Should this be considered a desirable 
property - which it certainly would be in a model where the 
immediate response to a step strain involved no plastic defor- 
mation - then this is a good reason to choose the midpoint 
treatment of the integral. 

5.1 Results for Delta-Step Approach 

For a single delta-step strain, the correlator obeys Eq. ( [T7] ) in- 
troduced previously. From this it is straightforward to estab- 
lish that the residual nonergodicity parameter, long after such 
a step strain, obeys 



f(i). 
no 



lim^«(0 



h{yo/2)mo[fo] 
l+/z(yo/2)mo[/o: 



From this follows the ultimate residual stress 



t(1)/ 



a^^^(^^oo)^a(oo)=yoVa(/4'^)^ 



(19) 



(20) 



which differs markedly in functional form from the ramp- step 
case, Eq. [15] Instead of saturating at large strains, as there, 
the residual stress now has a maximum at yo ~ Jc and falls to 
zero at large strain amplitudes, as was found p reviously in both 
ITT-MCT and schematic approaches'^^'^^^. (This result is 
also reminiscent of the Doi-Edwards theory and experiments 
on entangled polymer melts ^^.) The prediction is shown in 
Fig.[T]as a dashed blue line alongside data from the ramp-step 
approach. 

Through the fact that h decreases indefinitely for large 
strains, the delta-step model incorporates an extreme form of 
memory loss. This leads to plastic deformation in the time 
following the step strain that is very effective in reducing the 
stress. This contrasts with the ramp-step case (which, as ex- 
plained above, is closely related to the switchoff of steady 
shear): there, the memory loss during the application of the 
strain is self-limited due to the plastic flow already occuring 
in that time interval. 

Turning now to double step strain, the residual stress 
a(^) (oo) as a function of the two strain amplitudes is shown for 
reversing step strains in Fig. [5] where yo > but yi < in anal- 
ogy to the case discussed in Fig. [2] for the ramp model. The 
appearance of a maximum in the single- step response (dashed 




Fig. 5 Infinite-time stress plateau a(oo) after reversing double 
delta- step strains in the glass (F12 model, 8 = 10~^), as a function 
of the first and second step sizes, yo > and |yi | = — yi. Time 
between the strains: At = 10 (top graph) and A^ = 1 (bottom). 
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curve in Fig.[T]), reappearing here as the boundaries of the sur- 
face at Yo = or yi = 0, imposes a nonmonotonic variation 
of a*^^) (oo) also when both steps are finite. As noted in con- 
nection with Eq. ([18]), along the line yo = — yi, a(oo) ?^ is 
expected for small Al 

Along the line yo = — 2yi, a pronounced maximum occurs 
for small enough A^ that arises because of a cancellation of 
the yo/2 + yi factor in the damping function of Eq. ([18]). Inter- 
preting the one-sided 5-integral as unity would shift this max- 
imum to the line yo = — yi . No such maximum was found for 
the ramp calculations. It arises from a cancellation of strain- 
induced damping terms entering the memory functions and 
can hence be interpreted as an "echo" effect: two instanta- 
neous anelastic deformations can cancel if applied with a time 
lag Ar < To such that no intermediate Brownian motion de- 
stroys the reversibility. In the ramp model, on the contrary, 
plastic deformation is already present during the application 
of the strain steps, and these are not reversible. 

Just as in the ramp- step case (Fig. [2]), the plots in Fig. [5] 
show a strong departure from the form expected with BKZ- 
type constitutive equations. (Recall that in the latter, Eq. ([2]) 
evaluates to a(^)(oo) = 0^^:^^^ (oo), so that the graph in Fig. 5 
should be inflection- symmetric around the line yo = — yiT 
However the form of this departure is strongly different from 
that found in the ramp- step approach. For instance if ISt is 
large enough such that the response to the initial stress has 
decayed to its long-time plateau ay^^ (oo), the second strain su- 
perimposes a negative contribution such that the total resid- 
ual stress can become negative even if yo + yi > 0. But in- 
stead of a large second strain erasing all memory of the first 
(larger) step, the long-time asymptotes of the stresses essen- 



.(2) 



-^(l)r 



.«/ 



tially add: \im^t^ooO\^\j^{oo) = g\^\oo) +0^(00), i.e., the 
delta-step model keeps a persistent memory of the residual 
stress of the first step. 

For the case of exact strain reversal, results for the residual 
stress using the delta-step approach are shown alongside those 
for ramp- steps in Fig. [3] Apart from the much faster decay at 
large amplitudes (which is the counterpart of the very rapidly 
decreasing residual stress after a single step, discussed above) 
the two methods give broadly similar results. 



6 Conclusions 

For many strongly viscoelastic materials, such as polymers, 
the step-strain response is a cornerstone of nonlinear rheol- 
ogy. This holds both experimentally, where single and double 
step strains are major diagnostic tests, and in the development 
of theories. For instance, theories of BKZ type create an en- 
tire constitutive model solely from the nonlinear response to 
a single step strain, while double step strains have long been 
advocated specifically to probe any breakdown of BKZ-like 



superposition precepts (see, e.g. Ref.^^). Across the whole of 
polymer viscoelasticity theory, the assumption that step- strain 
responses can be measured experimentally, to sufficient accu- 
racy to inform these theories, goes largely unquestioned. 

In this paper we have presented a careful analysis of both 
single and double step strains within the schematic MCT ap- 
proach to the rheology of colloidal glasses. These are yield- 
stress materials, whose response to flow includes a strong con- 
tribution from plasticity: structural rearrangements are caused 
by straining the material. As discussed in the introduction, and 
amply confirmed by the calculations presented above, BKZ 
precepts can fail drastically in the description of such plastic 
materials. (This was also noted for a mesoscopic model of 
soft glasses in Ref. ^^ .) We have discussed this in terms of the 
residual stress a (00) remaining in a glass that has been subject 
to an exactly reversed double step strain. BKZ models pre- 
dicts this residual stress to vanish. Within MCT, a finite resid- 
ual stress remains: it is non- vanishing because plastic flow has 
taken place in the time between the strains, something not ac- 
counted for in the BKZ constitutive equation. This may have 
interesting consequences for other types of deformation like 
oscillatory shear. 

Perhaps more surprisingly, our work exposes the fact that 
in materials where plasticity is important, the very concept of 
a step strain is problematic. Specifically we have shown that 
within the mathematics of our schematic MCT model, the re- 
sponse to a discontinuous strain is different from the response 
to a finitely ramped one, even if the limit is then taken of the 
ramp rate tending to infinity. This is because in strong flows 
the relaxation time of the material is inversely proportional to 
strain rate. Accordingly, plastic deformation does not vanish 
in the fast ramp limit, whereas within our chosen equations it 
is (apparently) absent during the imposition of a mathemati- 
cally discontinuous step. 

A possible interpretation of the different representations of 
strain steps in schematic MCT is this: in the delta-function 
model, the strain in the system just after a step is primarily a 
large nonlinear elastic (anelastic) strain that then converts to 
plastic strain as memory loss kicks in. (Of course, this as- 
sumption is questionable whenever y^Jc, at least for hard- 
core potentials. It may be of relevance for soft particles.) The 
corresponding strain energy is relaxed via subsequent plastic 
deformations, as expressed through a strain-induced decay of 
correlation functions. The conversion of anelastic to deforma- 
tion to plastic rearrangements is particularly effective at large 
strains, so that most of the stresses relax, and a (00) is almost 
zero (dashed lines in Fig. [TJ. The strain-ramp model, on the 
other hand, assumes that plastic deformations are predominant 
already during the application of the strains. Consequently, 
the system is always fluidized in these time intervals, to the 
extent of completely erasing memory of any past flow history. 
While the delta- step model captures the fact that anelastic de- 
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formations can partially be transformed into plastic ones, but 
partially be recovered in a double- step strain setup, the ramp 
model captures that any strain much larger than y^ will lead to 
a full loss of reversibility. As a consequence, stresses can no 
longer be relaxed as effectively, and a saturation of the a(oo) 
curves is seen in Fig. [T] for Jo^Jc- 

Here, it should be noted, that the ramp model excludes some 
of the anelastic deformation mechanisms present in the full 
MCT. There the analog to Eq. ^ quite naturally leads to a 
vertex VaiJu') that itself depends on strain through a wave- 
vector advection mechanism^^^^. Furthermore, this vertex 
shows negative contributions for strains y ~ Yc, manifesting 
themselves in the appearance of a stress overshoot in startup 
flow^. Schematic models that include such anelastic contri- 
butions can be devised and will suppress the residual stress 
a(^)(oo) shown in Fig. [ijfor large yo, possibly reinstating a 
maximum around y^. 

Although we cannot be certain, the mathematical issue dis- 
cussed above might reflect a genuine physical one, of whether 
the bare Peclet number, yio "C 1 , is large or small during the 
rapid straining event that any experimental 'step' represents. 
(Here Xq is a non-glassy, microscopic relaxation time and y the 
shear rate.) In most 'step strain' experiments the bare Peclet 
number is unlikely to be large (see e.g. Ref.^); hence if this 
reasoning is correct, the results we have presented for finite 
ramp rates, rather than those for delta-function steps, are more 
pertinent. To confirm this point, it would be very interesting to 
see more experiments on step strain, particularly double step 
strain, in colloidal systems. As discussed above, the depen- 
dence of the residual stress a^^^ (oo) after a double step strain 
on the individual strains yo and yi can provide a sensitive test 
for the different scenarios. 

One has to keep in mind that real glassy materials will be 
subject to aging. In our discussion we implied that the initial 
unstrained configuration was equilibrated (i.e., representative 
of the Boltzmann ensemble, even if the dynamics is noner- 
godic in the glass). The usual rheology protocol (pre-shear to 
fluidize the system, followed by a certain waiting time ty^;) does 
not ensure this. We expect however, that for sufficiently long 
tw the schematic MCT results can at least qualitatively or even 
semi-quantitatively be compared to experiment, as was previ- 
ously found true for flow curves in steady state ^. Note that the 
imposition of a flow or step-strain displaces the system from 
Boltzmann equilibrium even if it was previously equilibrated 
as MCT assumes. The subsequent relaxation (e.g. between 
two steps in double step strain) may also have aging-like fea- 
tures, which MCT should capture. 
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A Appendix: Numerical Scheme 

For a strain history comprised of piecewise constant shear 
rates, one is led to split the integrals in Eqs. ([3j and ^ at 
the step discontinuities in y{t). Specifically, set y(^) = y^^^for 
tjn <t < tm+\ to obtain 

(21) 
for tn<t <tn^i. We have assumed y(^) = for t <to and have 
introduced functions G^^") {t^t') that are given by Eq. ^ with 
transient correlation functions (i;>''^^\t^t') defined for tm < t < 
tfn^i, tn<t' < tn+\. For thesc, the equations of motion read 

^htt' f htt''m[^^'^^\t-t'')]dt'4'^''\t\t')dt'' = 0. (22) 

Jtn 



Note that the correlation functions are transient ones, and as 
such are only affected by the strains applied between the times 
of their arguments. Hence, the (s^^'^'^\t,t') = ^^"^"^^t-t') are 
steady- state quantities (quiescent or steady- shear) that depend 
only on the time difference. Equations ([22]) are complemented 
by initial conditions obtained from the fact that correlation 
functions are continuous in both time arguments. 

We now specialize Eq. ([22]) to one or two strain ramps, 
y(0) = yq, y(2) = yj^ and all other Y^'") = 0. Then we are 
left with nine different correlation functions, three of which 
are steady-state ones: (^(^^^ = 4>*\ (^(^^^ = 0\ and ^(^^^ = 
(^(33) = ^^^ Figure [6] gives a schematic overview. It is straight- 
forward to specialize Eq. ^22) to the cases needed in evaluating 
a^^\t) and a^^\t), the single- and double-ramp response, 

^^^^ (0 =% f G!yl^^ (^ ^0 dt' for t > 8^ (23) 

G^'\t)=yol%f^%{t/)dt^ (24) 

+ ti / G\i (t,t') dt' for t>^t^ Iht. 

JAt^bt 



Here we have introduced interval lengths 5^ and At as indi- 
cated in Fig. [6] We can identify (|)(^^)(^,^') with the suitably 
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Fig. 6 Schematic view of the {t,t') plane indicating the various 
transient correlation functions appearing in the MCT equations for a 
double strain ramp. See text for details. 



shifted and extended ^^^^\t -At- bt,t' - At - 5t) evaluated 
for Yi instead of yo. The function (^''^^\t^t') defines the re- 
sponse to a single strain ramp, while (|)(^^) {t^t') includes cross- 
correlations between the two ramps. 

If one considers the limit 5^ -^ 0, the equations of motion 
for (|)(^^^ and (|)(^^) simplify as some of the integrals appearing 
in Eq. ([22]) vanish. Introducing x = t' /ht and ^{t^x) =^{t^t') 
as in the main text, we get for the single-ramp response 

[l+toar]C^(f,x) 

+ h{y(>{l-x))j"mMt-t")]d,4i''\t",x)=0, (25) 

a one-parameter family of functions for x G [0,1]. Similarly, 

(30) 



[i+toa,]CY,(^-^) 

rt 

+ h{yi+yo{l-x)) / mMt-t"M4^%{t",x)dt" = 

J At 

(26) 

There appears a cross-correlator in this equation, whose equa- 
tion (in the limit 5^ -^ 0) reads 

-hh{yif 1^' m[i^^^'\t,t')]d,.i^o{t' -t^)dt' 
+ /z(yi) / m[(^o{t-t'')]dt''^^^^\t'\t')dt'' = 0. (27) 

JAt 

This cross-correlator plays an important role in establishing 
the non-trivial dependence on At of the MCT double-ramp re- 



sults. It is a result of these and similar cross-correlations that 
the MCT constitutive equation is not of BKZ type. 

Deriving the equivalent set of equations for the case of 
delta- function steps, it is easy to see that again Eqs. ( [25] ) to 
( [27] ) appear, but evaluated only at one specific value of x; 
in the notation of the main text (^^^\t) = ^^^^\t,l/2) and 
(|)(^)(r) =^P^)(^,l/2). The cross-correlator (|)(^^)(^,^') appears 
identically in both models. 
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